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ABSTRACT 

We present simulated ./-band spectroscopy of red giants and supergiants with a 42 m European Extremely Large Telescope (E-ELT), 
using tools developed toward the EAGLE Phase A instrument study. The simulated spectra are used to demonstrate the validity of 
the 1.15-1. 22 fjm region to recover accurate stellar metallicities from Solar and metal-poor (one tenth Solar) spectral templates. From 
tests at spectral resolving powers of four and ten thousand, we require continuum signal-to-noise ratios in excess of 50 (per two- 
pixel resolution element) to recover the input metallicity to within 0. 1 dex. We highlight the potential of direct estimates of stellar 
metallicites (over the range -1 < [Fe/H] < 0) of red giants with the E-ELT, reaching out to distances of ~5 Mpc for stars near the tip 
of the red giant branch. The same simulations are also used to illustrate the potential for quantitative spectroscopy of red supergiants 
beyond the Local Volume to tens of Mpc. Calcium triplet observations in the /-band are also simulated to provide a comparison with 
contemporary techniques. Assuming the EAGLE instrument parameters and simulated performances from adaptive optics, the ./-band 
method is more sensitive in terms of recovering metallicity estimates for a given target. This appears very promising for ELT studies 
of red giants and supergiants, offering a direct metallicity tracer at a wavelength which is less afffected by extinction than shortward 
diagnostics and, via adaptive optics, with better image quality. 

Key words, instrumentation: adaptive optics - instrumentation: spectrographs - techniques: spectroscopic - galaxies: stellar content 
- stars: fundamental parameters 



1. Introduction 

Plans are well advanced for the next generation of optical and 
infrared (IR) ground-based telescopes, the extremely large tele- 
scopes (ELTs). Their science cases are broad, ranging from di- 
rect imaging of exoplanets, to studies of resolved stellar popula- 
tions in external galaxies, and spectroscopy of distant 'first light' 
galaxies at the highest redshifts (e.g. |Hook et al.|20 06). 

The ELTs are an increasingly global effort, with three 
projects under detailed study - the Giant Magellan Telescope 
(GMT), the Thirty Meter Telescope (TMT), and the European 
Extremely Large Telescope (E-ELT). With primary apertures in 
excess of 20 m, they will deliver a huge gain in our capabilities 
via a combination of unprecedented sensitivity and exquisite an- 
gular resolution. In parallel to the design of the observatories, 
significant effort has also been invested in studies of ELT instru- 
mentation (see |Jaffe et al.|2010t|Simard et al.|20l"0 ; Ramsay et al. 
[2QT0] >. 
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Over the past few years, deep imaging from ground-based 
telescopes and the Hubble Space Telescope (HST) has provided 
us with new and unique views of the outer regions of large galax- 
ies beyond the Milky Way for the first time, such as M31 (e.g. 
|Ferguson et"al~|2005] l and M33 (e.g. |Barker et al"~|20"07l >. From 
analysis of the resulting photometry and colour-magnitude di- 
agrams we can determine star-formation and assembly histo- 
ries for external galaxies, enabling tests of theoretical models 
of galaxy evolution (e.g. |Bullock & Jo hnston 2005), i.e. using 
resolved stellar populations as a tracer of the processes which 
have shaped the evolution of their host systems. 

The imaging from the Advanced Camera for Surveys (ACS) 
Nearby Galaxy Survey Treasury (ANGST; Dalcanton et al. 
2009 ) provides an excellent illustration of the diversity of galax- 
ies beyond the Local Group. Crucially, the Local Volume in- 
cludes a wide variety of morphological types - massive ellip- 
ticals, large metal-poor irregulars, lower-mass late-type spirals, 
interacting systems, dwarf starbursts - providing an excellent 
opportunity to quantify the effects of environment on galaxy evo- 
lution. 
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Photometric methods are powerful when applied to extra- 
galactic stellar populations but follow-up spectroscopy can aid 
our understanding significantly via precise chemical abundances 
and stellar kinematics. For example, results from the spec- 
troscopy of luminous blue supergiants in external galaxies ob- 
tained by the Araucaria project ( |Gieren et al.|2005 a). However, 
in pursuit of spectroscopy of evolved stellar populations, the 8- 
10m class telescopes are already near their limits beyond a few 
hundred kpc. For instance, observations with Keck-DEIMOS 
struggled to yield useful spectra below the tip of the red giant 
branch (TRGB) in M31 at I > 21.5 ( [Chapman et "aLl|2006> . If 
we aspire to spectroscopy of individual evolved stars in galaxies 
beyond the Local Group, we require the increased sensitivity of 
the ELTs, likely combined with some degree of correction from 
adaptive optics (AO) to mitigate the effects of crowding. 

Indeed, for stellar spectroscopy with the ELTs there will be 
a fine balance in sensitivity between the improved image quality 
from AO as one goes to longer wavelengths (where the wave- 
front errors become less significant compared to the observed 
wavelengths) versus the increased sky background. Moreover, 
work over the past decades has provided an excellent under- 
standing of many of the optical spectral lines available to us, 
enabling robust estimates of physical parameters and chemical 
abundances. To exploit the best performance (in terms of angular 
resolution) from the ELTs, we need to improve our knowledge 
of near-IR diagnostics, and also to pin-down the 'sweet spot' in 
terms of the gain in sensitivity from AO versus the background 
contribution. 

A recent study by Davies, Kudritzki & Figer (2010 here- 
after DKF10) suggested new /-band diagnostics, spanning 1.15- 
1 .22 fim (including lines from Mg, Si, Ti, and Fe), as a means 
to obtain stellar metallicities in extra-galactic red supergiants 
(RSGs). This wavelength region is relatively unexplored in this 
regard, only previously considered by Origlia et al. ([2004). 



Table 1. EAGLE baseline design. The patrol field is the instru- 
ment field-of-view within which integral-field units (IFUs) can 
be configured to observe individual targets/sub-fields. 



DKF10 noted that this method could be very powerful with 
new IR instruments under construction for 8-10 m class tele- 
scopes (such as VLT-KMOS and Keck-MOSFIRE), providing 
precise stellar abundances in galaxies out to a potential distance 
of ~10Mpc, complementing those from luminous blue super- 
giants (e.g. |Bresolin et al.|2001||Kudritzki et al.|2008) >. An even 
more compelling future prospect is in the context of ELT ob- 
servations, with DKF10 noting the potential of direct abundance 
estimates for individual RSGs in galaxies even beyond the Local 
Volume (subject to crowding). 

A second possible application of this spectral region is in 
observations of evolved red giant branch (RGB) stars - the long- 
lived descendants of much lower-mass stars. The marcs model 
atmospheres (Gustafsson et al. 2008) used by DKF10 are actu- 
ally of red giants (see their Section 2.2 for a discussion of the 
applicability of these models over a large range of stellar lumi- 
nosities). 

The potential of a direct abundance diagnostic for /-band 
ELT observations of both young (RSG) and old (RGB) stel- 
lar populations in external galaxies - with the benefit of bet- 
ter AO correction compared to shorter wavelengths - warrants 
further exploration. EAGLE is a conceptual design of an AO- 
corrected, multi-IFU, near-IR spectrograph, undertaken as one 



of the Phase A instrument studies for the E-ELT ( |Cuby et aL 
|2010| l. In this article we employ tools developed in the course of 
the EAGLE study as a proof-of-concept for quantitative /-band 
spectroscopy with the ELTs. 

We simulate EAGLE /-band observations with two objec- 
tives. First, to validate the technique of DKF10 for metal-poor 
spectral templates, and secondly, to explore the distances to 



Parameter 


Specification 


Patrol Field 


Eqv. 7' diameter 


IFU field-of-view 


1'.'65 x 1'.'65 


Multiplex (# of IFUs) 


20 


Spatial resolution 


30% EE in 75 mas (H band) 


Spectral resolving power (R) 


4000 & 10000 


Wavelength range 


0.8-2.5yum 



which we could obtain robust abundance estimates (for both 
RGB stars and RSGs) with ELT observations. To compare these 
results with diagnostics already in common use, we also con- 
sider simulations in the / band. Section|2]describes the tools and 
assumptions used in the simulations, with the analysis presented 
in Section|3] In Sections|4]and|5]we discuss the results in the con- 
text of the scientific potential of the E-ELT and other upcoming 
facilities. 

2. Simulations 

Three top-level science cases informed the requirements of the 
EAGLE design: the physics and evolution of high-redshift galax- 
ies, characterization of 'first light' galaxies at the highest red- 
shifts, and studies of resolved stellar populations in the Local 
Volume. The specifications of the baseline design are sum- 
marised in Table[T](with further discussion given by |Evans et al. 
2010b). 

To investigate the critical requirements for spatially-resolved 
spectroscopy of high-redshift galaxies, |Puech et al.| ( |2008| l de- 
veloped a web-based tool to generate simulated IFU datacubes. 



This was used by |Puech et al. (2008, 2010a) to investigate the 



required image quality, in terms of the level of AO correction, 
to recover spatially-resolved properties in high-redshift galax- 
ies (up to z ~ 6). Functional details of the websim simulator are 
given by Puech et al. (2010b). In brief, there is a web-interface to 
a core idl code. The user uploads input datacubes in the Flexible 
Image Transport System (FITS) format, which must incorporate 
specific keywords in the headers to provide information to the 
code, as well as sufficient spatial sampling (typically ten times 
greater than the size of the final IFU pixel) and the desired spec- 
tral sampling. The input cubes are then convolved with a model 
point-spread function (PSF), selected from a list of AO simula- 



tions (Section 2. 1 1 which are hard-wired into the code. By adopt- 



ing the relevant instrument/telescope parameters, |Puech et al. 



(2010a) provided an external check of the websim routines by 
reproducing datacubes to match those from SINFONI observa- 
tions of high-redshift galaxies. 

The parameters used in the websim simulations are sum- 
marised in Table|2] The current E-ELT design features a primary 
with an equivalent diameter of 42 m and with approximately 
9% obscuration by the secondary and supporting structures. The 
telescope throughput (>80%) was the expected temporal aver- 
age of the E-ELT design at the start of the EAGLE study, retained 
here as a worst-case estimate (compared to the >90% adopted by 
the Design Reference Mission studies, e.g., Puech et aL] 2010a). 

Inclusion of the sky background is discussed in detail by 
Puech et al. (2010a), in which a well-sampled model of the 
Mauna Kea background is scaled by measurements of the at- 
mospheric absorption at Paranal; the same approach is used 
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here. The E-ELT site was announced in early 2010 to be 
Cerro Armazones in northern Chile. Although at an altitude of 
3060 m, slightly higher than Paranal (the site of the Very Large 
Telescope, VLT), the adopted sky model is a good enough ap- 
proximation for our purposes. 

The websim routines add the sky spectrum (which includes 
the OH emission lines) to each spatial element of the datacube, 
then a noise contribution is added. In parallel, a separate sky 
spectrum is constructed using a different realisation of the noise 
(i.e. a different random number seed), to simulate an observed 
sky spectrum, which might be obtained by dedicated IFUs, on- 
off dithers, or using object-free spatial pixels in the same IFU. 
Each sky spectrum is then subtracted from the science spectrum 
for each spatial element, producing a sky-subtracted output cube. 
This technique is idealised slightly in the sense that there is little 
difference between the simulated sky and that in the science data, 
which could lead to residuals in the case of strong variations of 
the OH lines. However, OH residuals are present in the simulated 
spectra due to differences in the added noise to the science and 
sky components. 

In specifying the detector properties we adopted parameters 



Table 2. Parameters used in the websim simulations. 



Parameter 



from tests of Hawaii-2RG IR arrays by Finger et al. ( 2006 1: a 



dark current of O.Ole /pixel/s and a lowest read-out noise of 



2.3e /pixel (achieved using 'Fowler' sampling, see Fowler & 



Gatley 1990 1. Testing of the four Hawaii-2RG arrays destined 



for KMOS has yielded a read-noise of less than 10e /pixel from 
double-correlated sampling (Shar pies et al.|2010] >; Fowler sam- 
pling of the same arrays should approach the best-case results 
from Finger et al. with existing technology. Detector develop- 
ment is ongoing in parallel to plans for the ELTs (e.g. develop- 
ment of 4k x 4k arrays), such that detector cosmetics/properties 
will not limit ELT observations of the type that we describe here. 
Long detector integration times (1800s) were adopted, both for 
simplicity in the simulations and also to ensure that the observa- 
tions at R = 10000 are background limited (see Section [3.2.2| for 
further discussion). 

Our primary motivation is to explore the potential of the J- 
band spectral diagnostics, rather than a detailed trade-off study 
of spatial performance and required AO correction. In our sim- 
ulations we adopt the baseline EAGLE specifications for IFU 
slice-width (37.5 milliarcseconds), informed by the simulated 
observations of high-redshift galaxies Puech et al. (2008, 2010a). 
This leads to a (two spatial-pixel sampled) resolution element of 
75 mas on the sky. In terms of stellar observations this can be 
considered an intermediate level of AO correction, yielding ex- 
cellent image quality over a relatively large patrol field for effi- 
cient surveys. 

Note that many of the targets for absorption-line spec- 
troscopy with the ELTs are already known from deep HST imag- 
ing, but are beyond the reach of 8-10 m class facilities. Seeing- 
limited or Ground-Layer Adaptive Optics (GLAO) observations 
will not deliver sufficient contrast to match such imaging (e.g. 



Cunningham et al.|20 08 ), whereas the EAGLE sampling is com- 
mensurate with the HST pixel scale - i.e. if a target is resolved 
with HST, then EAGLE is well-matched to provide follow-up. 
The Near-IR Camera (NIRCam) for the James Webb Space 
Telescope (JWST) also has a similar pixel scale. 

In the densest regions in external galaxies or stellar clusters, 
finer angular resolution/greater contrast will be required (thus 
obtained over a smaller field). The spectroscopic methods dis- 
cussed here will be equally valid for instruments aimed at prob- 



Value(s) 



Diameter of primary (M] ) 42 m 

Effective central obscuration 9% 

Throughput (telescope) >80% 

Throughput (instrument, inc. detectors) 30% @0.850/jm 

35% @ 1.175 pm 

IFU slice width 37.5 mas 

Spectral resolving power 4000 & 10000 

Dark current 0.01 e~/pixel/s 

Read-out noise 2.3 e~/pixel 

Exposure time 20 x 1 800 s 



performance simulations would require tailored PSFs, beyond 
the scope of the current work. 



2.1. Simulated PSFs 

The EAGLE concept employs multi-object adaptive optics 
(MO AO) to provide significantly improved image quality for se- 
lected target fields within the focal plane (first advanced for the 
FALCON concept, e.g. |Assemat et al.|2007| >. The EAGLE base- 
line design uses an array of six laser guide stars (at a radius of 
3'75) and five natural guide stars (NGS) to map the atmospheric 
turbulence ( Rousset et al.|2010a I. The adaptive mirror in the E- 
ELT (M4) will be used primarily to correct the ground-layer con- 
tribution, while the higher-altitude turbulence will be corrected 
by deformable mirrors in each science channel. The baseline de- 
sign of EAGLE has an un-obstructed field-of-view with a diam- 
eter of 5', but the regions between the optics to monitor the laser 
guide stars are also accessible, giving a total field-of-view equiv- 
alent to a diameter of 7' . 

An extensive set of MOAO simulations were calculated as 
part of the EAGLE study ( |Rousset et aL|2010a|b i 



paying par- 
ticular attention to the location and magnitude of potential NGS 
(with R < 17) in the target fields. The PSFs used here were gen- 
erated using the ONER A wide-field AO analytical code ( |Neichel| 
et al.|200"8] ) with two real NGS configurations which, given the 
spatial distribution and magnitude of the available guide stars, 
are illustrative of the range of likely performance. These were 
initially selected for high-redshift simulations, so were taken 
for two example pointings within the area of the XMM-Newton 
Large-scale Structure Survey ( |Pierre et al. ||2004[), as shown in 



ing fin er scales such as HARMO NI for the E-ELT ( |Thatte et al 
2010) and IRIS for the TMT (Larkin et al. 2010), but detailed 



Figure [T] Five NGS were selected from those available in the 
'good' configuration (left-hand panel of Figure [T] which in- 
cludes one bright target with R = 10.5); only four NGS were 
available in the 'poor' configuration. As discussed by [Evans] 
|et al. I (2010b), the capability to use NGS beyond a 5' diameter 
leads to excellent sky coverage. 

In practise, the exact AO correction for each sub-field/target 
is a function of its location in the patrol field with respect to the 
NGS. Calculating tailored PSFs for each potential target is not 
practical, so we employ median PSFs (for both 'good' and 'poor' 
NGS configurations) derived from a grid of points across the 
EAGLE patrol field from the analytical code. Typical variations 
of the ensquared energy (EE) across the patrol field are of order 
±10% - a key part of the planning for such observations will be 
to optimise the selection of targets in terms of the expected AO 
performance. 

Assuming a 'Paranal-like' set of conditions, we considered 
two sets of seeing in the simulations (with the same vertical tur- 
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Science Field: XMM-LSS, sub-field 10 



Science Field: XMM-LSS, sub-field 1 
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Fig. 1. Configurations of potential natural guide stars (NGS) used in the simulated PSFs (labelled by their /?-band magnitudes). The 
overlaid circles have diameters of 5', 7.3', and 10'. Left: 'good' NGS configuration, with multiple stars available. Right: 'poor' NGS 
configuration, with only four stars within the intermediate circle. 



bulence profile): 0765 at A — 0.5/im at zenith (the mean VLT see- 
ing; Sarazin et al.| 2008 ) and 0790 at a zenith distance (ZD) of 
35°, illustrative of the performance that one might expect from 
execution of a 'Large Programme' -like survey. The EE in 75 mas 
for each PSF is summarised in Table [3] 

These performances are knowingly optimistic in the respect 
that fitting, tomography and noise errors were included, but ad- 
ditional error terms will contribute to the final performance. 
|Rousset et al.| ( |2010b[ ) evaluated the additional wavefront error 
from these aberrations, which account for the limitations of the 
tomographic model, calibration, and operating in open loop. The 
EE degradation can be approximated by a factor of e~°~^~, where 
crip 1 is the phase variance (in radian 2 ) of the wavefront error at 
the observed wavelength. This leads to a scaling of the EE in the 
simulated PSFs by a factor of 0.69 in H (calculated at 1.65 fim), 
0.52 in / (1.25 fan), and 0.25 in / (0.85 /an). Although the mor- 
phology of the PSF will also change slightly, in calculation of 
the resulting signal-to-noise (S/N) ratios in the simulated spec- 
tra (next section) and interpretation of our results (Section |4| 
we consider only the more important equivalent loss in sensi- 
tivity. As noted before, the current study is primarily intended 
as a proof-of-concept study for the /-band diagnostics - the ex- 
act AO performances will continue to evolve with the telescope 
and instrument design, and experience from on-sky demonstra- 
tors such as CANARY on the William Herschel Telescope (e.g. 
|Gendronetal.|2010| >. 



2.2. Template spectra 

Our simulations adopt two spectral templates from the grid of 
low-temperature stellar atmospheres calculated with the marcs 
code (Gustafs son et al.| 2008). These have an effective temper- 
ature (T e ff) of 3800 K (corresponding to an equivalent spectral 
type of approximately M0), log g = 0.0, and a microturbulence 
(£) of 2 km s . Models for two metallicities (Z) were used, with 
log[Z] =0.0 and —1.0 (normalised to Solar). The latter model 
(i.e. 10% Solar) serves as a test of our methods in the low 
metallicity regime, comparable to, e.g., results for some of the 
young supergiants observed in the WLM dwarf irregular galaxy 



Table 3. Ensquared energy (EE) in 75 milliarcseconds in the 
simulated MOAO PSFs used in the /- and /-band simulations 
presented here. Note that these do not include some of the po- 
tential error terms (see Section |2~T} , and are scaled accordingly 
in the discussion. 



Band 


A c [pm] 


Seeing 


ZD 


NGS 


EE [75 mas] 


/ 


0.85 


a/65 


0° 


Good 


43% 


/ 


0.85 


«' 65 


0° 


Poor 


37% 


/ 


0.85 


a.' 90 


35° 


Good 


25% 


/ 


0.85 


a/ 90 


35° 


Poor 


19% 


J 


1.25 


a/ 65 


0° 


Good 


64% 


J 


1.25 


a/ 65 


0° 


Poor 


58% 


J 


1.25 


a.' 90 


35° 


Good 


48% 


.1 


1.25 


a/ 90 


35° 


Poor 


40% 



( Urbaneja et al.||2008 l. Further details of the marcs models, in 
the context of the current study, are given by DKF10. 

At Solar metallicity we also calculated a set of /-band sim- 
ulations to inform discussion of the /-band results in Section [4] 
The /-band includes the calcium triplet (CaT, spanning 0.85- 
0.87/mi), which has become an increasingly ubiquitous diagnos- 
tic of stellar metallicities and radial velocities in nearby galaxies 
(e.g. Tolstoy et al. 2004). The CaT is an empirical calibration 
of metallicity (which employs the equivalent width of the Ca n 
lines) rather than a direct diagnostic. However, careful calibra- 
tion yields a robust relationship over a large range of metallicity 
(-2.5 < [Fe/H] < -0.2; |Cole et al.|2004) . Our primary interest in 
the CaT is simply to obtain an estimate of the continuum sensi- 
tivities, so only one template (log[Z] = 0.0) is required. 

The templates were convolved by the relevant spectral re- 
solving power, and then binned such that each resolution ele- 
ment is sampled by two pixels. The templates were restricted 
to spectral windows of AAl. 15-1.22 /an (/-band) and /1/10.840- 
0.875 //m (/-band), with artificial continuum regions appended 
(with an intensity set to unity) to provide straightforward esti- 
mates of the S/N in the simulated spectra. 

The websim routines operate on fluxes (F) in physical units 
(either ergs/s/cm 2 /Hz or ergs/s/cm 2 /A), so the input model tem- 
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Table 4. Summary of continuum signal-to-noise (S/N) obtained per two-pixel resolution element for simulated /-band spectroscopy 
(fexp = 10 hrs, R = 4000 & 10000). The quoted magnitudes are after scaling by the estimated uncertainties in the simulated PSFs (see 
Section 2.1 for further discussion and for details on the configurations of natural guide stars, 'NGS'). 
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44+3 


35 + 2 


55+4 


49+2 


23.25 


0.0 
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30+ 2 


23 + 2 


38+2 


36+ 2 


23.75 


0.0 


35 ± 


4 


27+3 


45+6 


40+ 3 












-1.0 


32 ± 


5 


27+4 


40+ 4 


39+ 5 











plates were scaled to the desired observed magnitudes (M, in the 
Vega system), using: 



F v = 



"M„=0 



1Q(M,/2.5) 



(1) 



where F, 



Campins et al. 



0) : 



1985 



2550 Jy and Frj=o) = 1600 Jy (|Bessell||1979 
respectively^] 



The simulated IFU observations are of an isolated point 
source centred in the central (37.5x37.5 mas) spatial pixel, 
which is then extracted to yield the object spectrum for each tar- 
get magnitude/set of conditions. The image quality of the sim- 
ulated PSFs is sufficiently good that, in this idealised case of a 
centrally-located point source, extraction and merging with adja- 
cent spatial elements does not improve the S/N of the final spec- 
trum; in real EAGLE observations, more sophisticated extraction 
routines (e.g., PSF-fitting) would likely be required. 



Table 5. Summary of continuum signal-to-noise (S/N) obtained 
per two-pixel resolution element for simulated /-band spec- 
troscopy (f e xp = 10hrs, /? = 10000). As in Table [4j the quoted 
magnitudes are after scaling by the estimated uncertainty in the 
simulated PSFs (see Section [2~T| . 





Seeing = 079 @ ZD=35° 


Seein 


I = 0765 @ ZD=0° 


1\EGA 


NGS good 


NGS poor 


NGS 


good 


NGS poor 


20.1 


121 + 6 


86 + 4 


177 


±9 


150 + 7 


20.6 


89 + 4 


64 + 2 


133 


±8 


111+7 


21.1 


64 + 3 


45 + 2 


101 


±2 


83 + 3 


21.6 


46 + 2 


30 + 2 


72 


±3 


60 + 2 


22.1 


31 + 2 


21 + 1 


51 


±2 


42 + 2 


22.6 


20+1 


14+ 1 


36 


± 1 


29+1 


23.1 


14+ 1 


9+1 


24 


± 1 


19+1 


23.6 


9+1 


6 + 


16 


± 1 


13 + 1 



2.3. Simulation results 



Each set of simulations was run ten times to provide a mean- 
ingful estimate of the dispersions in the derived metallicities and 
continuum S/N. Example simulated /-band spectra are shown in 
Figures [2] and [3 with notable diagnostic lines of Mg i, Si i, Ti i 
and Fe i labelled. Example /-band spectra are shown in Figure[4] 
The means and standard deviations of the continuum S/N 
from the /- and /-band runs are summarised in Tables [4] and 
[51 respectively; the results are also shown in Figures |5j [6] and 
Pj Following the discussion in Section 2.1 regarding the known 
limitations of the AO simulations, the quoted magnitudes take 
into account the factor of two in intensity in the /-band (equiva- 
lent to A/ = -0.75) and the factor of four in the /-band (equiva- 
lent to Al - -1.5), i.e. the effective magnitude of each simulation 
is brighter than the input magnitude. All quoted S/N results are 
for spectra extracted from the central spatial pixel of each cube, 



1 The passband of Cousins / is somewhat different to that of the 
F8 14W HST filters, but the typi cal magnitude offset is small: / - F8 14W 
< 0.05 m ( (Holtzman et al.|1995) . 



per two spectral pixels (set to be the sampling of the resolution 
element in each A and R combination). 



3. Analysis 

3.1. J -band simulations 

Our goal in the first part of the analysis is to determine physical 
properties from the simulated observations and compare them 
with those of the input model, i.e., how successfully can we re- 
cover accurate metallicities at both log[Z] =0.0 and -1.0, and 
what S/N is required? 

We follow similar methods to those described by DKF10. 
Each simulated spectrum is compared with template spec- 
tra from the marcs grid, the x 1 is calculated, and the loca- 
tion in the grid of the x 1 minimum is found. The parameter 
space searched includes 3200 < T e g < 4000 K (in steps of 200 K), 
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Fig. 2. Examples of simulated fl = 4000 spectra for 7 = 21.25, 22.25, and 23.25 (with seeing of 079, ZD = 35°, and the 'good' 
NGS configuration). Upper panel: Solar metallicity simulations (log[Z] =0.0); lower panel: metal-poor simulations (log[Z] = -1.0). 
Identified lines, from left-to-right by species, are: Ti i /LU.1896, 1.1953yum; Si i /t/U.1988, 1.1995, 1.2035, 1.2107 /mi; Mg i 
/Ltl.1831, 1.2087//m;Fe 1/1/11.1597, 1.1611, 1.1641, 1.1693, 1.1783, 1.1801, 1.1887, 1.1976//m. Residuals from the sky OH lines 
can be seen in places, e.g. the emission 'spikes' in the simulated metal-poor spectrum for J = 23.25. 



- 1 .5 < log[Z] < + 1 .0 (in steps of 0.25 dexf] 2 < <f < 5 km s" 1 (in 
steps of 1 km s ^ and log g = 0.0 and +1.0. 

Two refinements to the methods from DKF10 were neces- 
sary for the analysis presented here. This concerns the location 
of the continuum level in the simulated spectra and the interpo- 
lation between model grid points to refine the metallicity deter- 
minations, both of which are now discussed. 

3.1.1. Continuum placement 

To obtain accurate metallicities from the simulated spectra it is 
crucial that the continuum level is determined correctly. Placing 
the continuum at the wrong level results in an incorrect measure- 
ment of the normalised line strengths, thus an erroneous metal- 
licity. DKF10 found the continuum level of their input spectra 



2 The grid does not contain models for log[Z] = -1.25, so spectra for 
this metallicity are interpolated from neighbouring models. 

3 The spectra at £ = 3 and 4 km s~' are interpolated (see DKF10). 



by ranking the spectral pixels in order of their (normalised) flux, 
and then finding the median value of those pixels with the great- 
est flux. In a spectrum comprised of continuum and absorption 
lines, the pixels with the greatest flux will have a value equiva- 
lent to the continuum level plus the noise. At high S/N this is a 
very good approximation of the continuum level, but at low S/N 
this can lead to the estimated continuum level being too high, re- 
sulting in a derived metallicity that is systematically above that 
of the input spectrum. 

Here we adopt an improved method of continuum determina- 
tion, assuming that the noise is Gaussian (i.e. requiring good sky 
subtraction). A spectrum of pure Gaussian noise, when ranked 
in order of increasing pixel value, will produce a trend which is 
approximately linear in the middle-ranking pixels (see Figure 8}. 
This linear regime is selected to contain the values that are within 
~lcr of the mean. An absorption-line spectrum (i.e. discrete fea- 
tures combined with Gaussian noise) will also have this linear 
regime (see Figure|9]l. By identifying and fitting these pixels, we 
can find those within some tolerance level of the fit. In princi- 
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Fig. 3. Examples of simulated R = 10000 spectra for 7 = 21.25, 22.25, and 23.25 (with seeing of 079, ZD = 35°, and the 'good' 
NGS configuration). Upper panel: Solar metallicity simulations (log[Z] =0.0); lower panel: metal-poor simulations (log[Z] = -1.0). 
Identified lines are the same as those in Figure|2] 



/ = 20.6 : 




8450 8500 8550 8600 

Wavelength (Angstroms) 



8650 



8700 



Fig. 4. Examples of simulated CaT spectra (R = 10000, log[Z] = 0.0) for / = 20.6, 21.6, and 22.6, (with seeing of 079 and the 'good' 
NGS configuration). The three Ca n lines identified are AA8498, 8542, 8662 A. 
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pie we can then select those within ~l<x of the mean continuum 
level, and use them to fit the continuum. Figures [8] and ^illus- 
trate this procedure for a spectrum of pure Gaussian noise and 
for a template absorption-line spectrum with noise added (both 
withS/N = 50). 



3.1.2. Refinement of the metallicity measurement 

The marcs grid provides models at intervals of 0.25 dex between 
log[Z] = —1.5 and +1.0 (excepting -1.25). Determining a finer 
estimate of the stellar metallicity requires that we interpolate 
within the model grid. DKF10 achieved this by taking a x 2 - 
weighted mean of the best-fitting points in the grid. However, 
if the x 1 trend with metallicity is asymmetric about the x 2 min- 
imum (i.e. the best-fitting value), the 'refined' metallicity will 
always be systematically offset from the best-fitting value in 
the direction of the shallowest x 2 gradient. For example, if the 
X 2 minimum in the grid is at log[Z] =0.0, but the x 2 diverges 
more rapidly at negative log[Z] than at positive values, the^ 2 - 
weighted mean metallicity will always be above 0.0. Indeed, in 
the course of the current study we found that determining stellar 
metallicities in this way produced 'refined' estimates that were 
systematically larger than the input values (a consequence of the 
morphology of the x 2 -minimum). 

To eliminate these systematic errors we introduced a new 
method to refine the metallicity determination beyond the res- 
olution of the model grid. After finding the location of the x 2 
minimum in the grid, we extract a two-dimensional plane in T e ff- 
log[Z] space at equal log g and £ values. We then resample this 
grid onto a grid five times finer using bi-cubic spline interpola- 
tion, and take the location of the^- 2 minimum of the resampled 
grid as the refined metallicity and stellar temperature. 



3.2. Results 

Each simulated spectrum was run through the atmospheric anal- 
ysis code to determine r e jf, log[Z], £ and log g. The mean metal- 
licities recovered (and their standard deviations) from the ten 
runs of each simulation (for input metallicity, R, NGS configu 
ration, seeing, and effective magnitude) are shown in Figure 10 



3.2.1. # = 4000 

The mean metalliticites recovered from the R = 4000 simulations 
have a small scatter around the input value (<0. 1 dex) down to 
/ = 22.75 (see Figure 10 1. This is the point at which the errors 



start to become larger, although the dispersion is only ±0.1- 
0.15 dex for / = 23.25 with a/65 seeing. 

These results provide a good test of the DKF10 methods for 
metal poor stars, demonstrating reliable metallicity estimates at 
both metallicities above a given S/N threshold. Of course, S/N is 
not the only effect which determines the accuracy of the metal- 
licity estimate. Other effects will likely play a role, such as the 
robustness of the atomic data, any degeneracy between temper- 
ature, surface gravity, and microturbulence, and which part of 
the curve-of-growth the selected lines are on (i.e. the effects of 
saturation on the observed equivalent width and its relation to 
abundance). Nevertheless, comparing these results with the S/N 
values in Table|4]suggests a minimum S/N of 50 to 60 is required 
to recover Z to within 0. 1 dex. 



3.2.2. R = 10000 

Analysis of the R = 10000 simulations also recovers good metal- 
licity determinations (AZ + 0. 1 dex) for both templates. The min- 
imum S/N threshold is marginally lower than for the /? = 4000 
spectra, with S/N > 45 required from comparison with the results 
in Table g] 

Given the greater spectral dispersion, the faint magnitude 
limit is (at least) half a magnitude brighter than for the R = 4000 
simulations. Thus, use of the lower resolving power is (as to be 
expected) the most sensitive option if just a global estimate of 
metallicity is required for a given target. However, if precise ra- 
dial velocities are also required, to map the dynamics of a given 
region of the host galaxy, then observations at R= 10000 pro- 
vide improved velocity resolution (30 km s , as compared to 
75 km s- 1 at fl = 4000). 

To investigate the noise properties in the higher-resolution 
mode, the intermediate files from the websim were analysed. 
The sky background continuum counts (~60 photons/pixel for 
integrations of 1800s) give a shot noise that is over three times 
greater than the adopted detector read-noise. The noise on the 
dark current is actually the more critical factor in this regard, 
but note that our adopted dark value is somewhat conservative 
(e.g., a factor of two greater than that measured for the arrays 
tested for JWST by Figer et al. 2004). Shorter integrations of 
900s remain background limited, although minimising the detec- 
tor read-noise in this instance becomes an even more important 
requirement. 



3.3. I-band simulations 

Understanding the relative performances for /- and /-band ELT 
spectroscopy will be valuable when planning future observa- 
tions. Exploratory simulations of the CaT were given by [Evans] 
etal.| (2010a) to investigate the potential of EAGLE spectroscopy 
in this region. To inform discussion of the new /-band results, 
we revisited simulations of the CaT, now adopting the same 
marcs templates and seeing/NGS configurations as for the J- 
band simulations. The adopted instrument throughput is slightly 
lower at 0.85 jum (30%) compared to 1.175 /urn (35%) in the fi- 
nal EAGLE design. As before, these values are not adopted to be 
definitive, but are representative of a near-IR instrument solution 
in which the primary performance requirements are at longer 
wavelengths. Table[5]summarises the mean continuum S/N (and 
standard deviations), per two-pixel resolution element, obtained 
from sets of ten simulated spectra calculated at half-magnitude 
intervals^] 

As an external check on the simulated /-band performances 
we turn to the ESO E-ELT spectroscopic exposure time cal- 
culator (ETC^j Input parameters were adopted to match those 
in the EAGLE simulations as closely as possible: a marcs 
model atmosphere with r e ff = 4000, point-source target, Paranal 
site, 42 m primary, seeing of 078, zenith distance=30°, laser- 
tomography/multi-conjugate AO PSF, R= 10000, a 20 mas ra- 
dius of the S/N reference area (i.e. best match to our slice- 
width of 37 mas), and individual integration times of 1800 s. For 
S/N > 10 (per two-pixel resolution element), the ETC requires 
10 hrs of total integration time for / = 24.85. Note that the simu- 



4 The simulations were originally calculated adopting the 
/-band zero-point from the //Sr-NICMOS conversion tools: 
http://www.stsci.edu/hst/nicmos/tools/conversion_help.html (2250 Jy at 
0.9/iin). Translated to the Bessell (1979) system this results in a shift 
of 0.1 m fainter, hence the sampling in Table|5] 

5 http://www.eso.org/observing/etc/ 
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lated PSFs incorporated in the ESO ETC have some of the same 
limitations to those generated for EAGLE, taking this into ac- 
count leads to an effective /-band magnitude of ~23.5, i.e., of 
comparable sensitivity to the results in Table [5] Given the dif- 
ferences in the AO architectures, the exact modelling tools (and 
assumptions therein), and the methods used to calculate the S/N, 
the ETC comparison provides an independent check that our 
sensitivity limit (for a given S/N) is within the expected range. 



4. Discussion: RGB stars 

We now consider the simulation results in the context of the 
(unreddened) distances to which we can recover accurate metal- 
licities from /- and /-band spectroscopy with the E-ELT. For 
reference, absolute /-band magnitudes (Mj) for RGB stars span 
-1 <M 7 <-3, with typical colours of (/ - J) ~ 0.7-0.8 (e.g. 
Girardi et aLpOO^ |Marigo et al.|2008) . At the tip of the RGB 
absolute magnitudes extend up to Mj = -4. 



4.1. Potential for E-ELT & EAGLE 

We require S/N > 50 for accurate metallicities from the /-band. 
At R = 4000, the results in Table [4] provide sufficient S/N down 
to J -22.15 to 23.25, depending on the NGS configuration and 
seeing. Taking the bright end of the RGB (adopting Mj = -3.75 
from the colours above), this corresponds to an unreddened dis- 
tance modulus of 27, or a distance of ~2.5 Mpc. 

Thus, the DKF10 methods could provide a direct metallic - 
ity indicator in RGB stars out to a volume which encompasses 
a diverse range of galaxies. For instance, as we move beyond 
the Local Group, there are 32 known galaxies with distances in 



the range 1.0 < d < 2.65 Mpc ( |Karachentsev et al.||2004[ >. Most 
notably these include NGC 55 and NGC 300, the two spirals in 
the closer part of the Sculptor 'Group' at 1.9 Mpc (Gieren et al. 
|2005b| [20581 |Pietrzynski eTdl|2006) . NGC 55 is seen to have 
'LMC-like' abundances from spectroscopy of luminous super- 
giants ( |Castro et al.|2008| & in prep) and NGC 300 has a radial 
oxygen gradient that ranges from 0. 1 dex below Solar to near 
LMC levels at a radius of ~5 kpc (Bresolin et al. 2009 and refer- 
ences therein). This range of distances also includes lower-mass 
galaxies such as Sextans A, NGC 3109, and GR8, which are all 
depleted in metals (as traced by oxygen) by approximately 1 dex 
UKaufer et al.|2004[ |Evans et al.|2007) |Van Zee et aI1[20^6l re- 
spectively). Of course, some galaxies within 2.5 Mpc, in addi- 
tion to the faint dwarfs in the Local Group, will have even lower 
metallicities, providing motivation to further explore these diag- 
nostics in the very metal-poor regime. 

Scaling these results for TRGB stars (Mj — 5), good metal- 
licities should be recoverable out to distances of 4 to 5 Mpc, 
depending on the exact AO correction. This opens-up an 
even wider range of targets, encompassing over 200 galaxies 
( Karachentsev et al. 2004 ), including those in the Centaurus A 



Group (e.g. Karachentsev et al. E002). Direct metallicities for 



TRGB stars will play an important role in calibration of the dis- 
tance scale in the local universe, reducing the uncertainties asso- 
ciated with photometric distance determinations from the TRGB 



method (see discussion by Kudritzki 


2010), which can lead to 


systematic uncertatinies of ~5% (e.g. 


Mager et al.|2008 1. 



Table [5] To determine the limiting magnitude for metallicities 
from our /-band simulations, we adopt a requirement of S/N > 20 
(per two-pixel resolution element). Battag lla et al.| ( |2~008| l con- 
sider S/N > 20 (per A) as the threshold to obtain metallicities to 
within 0. 1 dex from CaT spectroscopy at R = 6500 (from VLT- 
FLAMES). At the centre of our /-band template 1 A is equiv- 
alent to 2.33 pixels, so the S/N results in Table [5] are slightly 
conservative in this regard, but they are also at greater spectral 
resolving power. However, the dominant factor here is the lower 
AO performance in the /-band than at longer wavelengths. The 
EE in the simulated PSFs is already considerably lower than in 
the /-band (Table[3]l, but the additional error terms are also larger 
(i.e. a factor of four compared to a factor of two in the /-band). 

Depending on the seeing conditions and NGS configuration, 
the /-band simulations yield S/N > 20 for 22.1 </<23.1, com- 
parable to the /-band limit discussed in the previous section 
(S/N > 50 at /? = 4000). Given that RGB/TRGB stars are mod- 
erately red (/ - / = 0.6-1.0), the DKF10 technique is probing 
deeper in terms of absolute luminosity of a given star. 

The /-band technique exploits the improved AO correction 
at longer wavelengths. It requires greater S/N to recover a good 
estimate of stellar metallicity but, given the final EAGLE instru- 
ment throughputs and current expectations of AO performance 
shortwards of 1 yizm, outperforms the CaT observations for a 
fixed exposure time of a given RGB star. In addition to a direct 
metallicity measurement (rather than the CaT calibration, albeit 
thought to be well understood), /-band observations will have 
the advantage of reduced extinction compared to observations 
at shorter wavelengths. Although we have limited ourselves to 
consideration of extra-galactic sources, working at longer wave- 
lengths would be attractive for observations in, for example, ob- 
scured Galactic clusters. 



4.3. Potential for TMT-IRMS 

Our simulations can also be considered in the context of the 
Infra-Red Multi-Slit Spectro meter (IRMS) being developed for 
the TMT ( |Simardetal.|2010| ), which is a clone of the MOSFIRE 
instrument for the Keck I telescope. IRMS will be located behind 
the TMT multi-conjugate AO system (NFIRAOS), taking ad- 
vantage of excellent image quality across a 2' field. Its planned 
wavelength coverage is 0.9-2.5 fim and, depending on the slit- 
widths, the delivered resolving powers will be R = 3000 to 5000. 

The minimum (Nyquist sampled) IRMS slit-width will be 
160 mas, i.e. coarser spatial sampling and with a greater sky 
contribution for a point-source target than from EAGLE. The 
expected /-band encircled energies for a slit- width of 160 mas 
range from approximately 30 to 45 %, depending on the posi- 
tion within the 2' field (Dr. Brent Ellerbroek, private communi- 
cation). These simulations include some of the additional error 
terms factored into the discussion of our MOAO PSFs, i.e., when 
the larger IRMS slit-widths are taken into account, they provide 
a comparable range of image quality. 

Assuming that the sensitivity varies as a factor of primary 
aperture (combined with a greater background flux from a larger 
effective slit), this corresponds to observations of RGB stars out 
to distances of ~1.5 Mpc and stars at the TRGB out to ~3 Mpc 
(depending on limitations of the NFIRAOS performance esti- 

Tailored performance 



4.2. Comparison with I-band performances 

To compare the /-band results with the commonly used CaT 
diagnostic, we now consider the /-band performances given in 



mates cf. those discussed in Section 2.1 
simulations are warranted for stellar spectroscopy specific to the 
TMT design, instrumentation and site (selected to be Mauna 
Kea in Hawaii). Nevertheless, we highlight the potential of the 
DKF10 methods for TMT spectroscopy of stellar populations 
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beyond the Local Group. For example, in NGC 1569, a north- 
ern starburst galaxy at 1.95 Mpc ( Karachentsev et aL|2 004 ) with 
a metal abundance that appears to be intermediate to those of 
the LMC and SMC ( |Devost et al.|1997,[Kobulnicky & Skillman| 
1997, cf., for exampleT TTrundle et al.|2007) . The TMT will also 
be well placed to explore the wealth of faint targets in northern 
hemisphere galaxies that are closer to home, e.g., M31 and M33. 



5. Discussion: Red supergiants 

The work by DKF10 was originally motivated by spectroscopy 
of RSGs. With the benefit of AO correction and their intrinsi- 
cally red colours, ELT observations of RSGs will provide our 
most distant quantitative tracers of individual 'normal' stars. 
Extending the earlier discussion, and with absolute magnitudes 
of-7<M/<-ll, 7 = 23.25 (atR = 4000) corresponds to an im- 
pressive 70 Mpc for the most luminous RSGs. The potential of 
such observations requires investigation in greater detail. In par- 
ticular, the impact of crowding (and the required AO correction) 
will need to be assessed, for instance, using the HST observa- 
tions of galaxies in the Virgo and Fornax Clusters (Cote et al. 
|2004l|Jordan et al.|2007] >. 



5.1. Potential for VLT-KMOS 

More iminently, these results are relevant for KMOS (under con- 
struction for the VLT), which will have 24 IFUs, each 278 x 278 
on the sky (with 072/slice), providing R = 3400 in the YJ bands 
(Shar pies et al.|20T0] l. The ETC for the natural-seeing mode of 
SINFONI provides the most appropriate comparison available 
at present - the limiting magnitude for S/N ~ 50, in a total expo- 
sure time of lOhrs, is / = 18.75. As noted by DKF10, this will 
enable us to obtain surveys of RSGs in galaxies well beyond the 
Local Group, and for the brightest targets at distances of up to 
~ 10 Mpc. This will provide the first opportunity to obtain large 
samples of RSGs to investigate large-scale metallicity gradients 
and structure in external galaxies. 



5.2. Potential for JWST -NIRSpec 

These methods will also have applications with the JWST Near- 
IR Spectrometer (NIRSpec), which is a multi-object spectro- 
graph capable of observing over 100 sources over a field of 
approximately 3.5x3.5 arcmin. NIRSpec will have three se- 
lectable spectral resolving powers (R = 100, 1000 & 2700), over 
a range of 1 -5 //m. 

The expected point-source NIRSpec sensitivit^for S/N = 10 
(two-pixel resolution element) with the R = 2700 grating, at 
1.15//m, is 950.7 nJy (equating to /vega = 23.1) in a total in- 
tegration of 100 ks (27.8 hr). This suggests a limiting magnitude 
of / = 19.6 (for S/N ~ 50) which, combined with the large mul- 
tiplex of the mult-slit array, could obtain large samples of RSGs 
in galaxies out to the Virgo Cluster. 



6. Summary 

We have simulated /-band EAGLE/E-ELT observations of an 
M0 giant (and, by extrapolation of the model atmospheres, an 
M0 supergiant). We extended the work of DKF10 for the 1.15- 
1 .22 fim diagnostic region, finding: 



http://www.stsci.edu/jwst/instruments/nirspec/ 



- Robust stellar metallicities can be obtained from analysis of 
/-band spectra of metal-poor stars (log[Z] = —1.0). 

- Accurate metallicities (±0.1 dex, at log[Z] = -1.0 and 0.0) 
require: S/N > 55 at R = 4000, and S/N > 45 at R = 10000. 

- The Phase A design of EAGLE on a 42 m E-ELT has the 
potential of direct metallicity estimates for RGB stars out to 
distances of ~2.5 Mpc, and for RSGs well beyond the Virgo 
and Fornax Clusters (subject to crowding). 

- For our instrument assumptions and simulations of rela- 
tive AO performance, the /-band diagnostics are slightly 
more sensitive than the CaT (A mag ~0.5-1.0) for recover- 
ing metallicities of a given MO-type star. The /-band method 
also has the advantage of reduced extinction compared to the 
CaT region. 

The first two of these points are not specific to EAGLE (nor 
dependent on the simulated PSFs used here), serving as a proof- 
of-concept for broader application of the methods advanced by 
DKF10. This part of the /-band is relatively free of OH emission 
lines compared to other parts of the near-IR, enhancing its po- 
tential for quantitative spectroscopy of extragalactic stellar pop- 
ulations with ELTs. This diagnostic range also has great poten- 
tial for other upcoming facilities such as VLT-KMOS, JWST- 
NIRSpec, and TMT-IRMS. 

One of the prime drivers for inclusion of the /-band in the 
design of near-IR spectrographs for ELTs is to study extragalac- 
tic stellar abundances via the CaT. If comparable or better re- 
sults can be obtained from the /-band, inclusion of the /-band 
would become less critical, leading to simpler and potentially 
cheaper instrument designs (e.g. fewer gratings), and avoiding 
compromises which may detract from the potential performance 
at longer wavelengths. 

As noted by DKF10, an important next step is quantitative 
analysis of RSGs in the Magellanic Clouds to refine these tech- 
niques over a range of metallicities, and to also extend this work 
to RGB stars in the Clouds. These methods should also be ex- 
plored at yet lower metallicities (relevant to the study of metal- 
poor galaxy halos), both via analysis of metal-poor observations, 
and by use of a larger grid of theoretical calculations. 

Diagnostics in the H- and /f-bands (in which the AO cor- 
rection will be even better) also warrant consideration in this 
context. For instance, use of /if -band equivalent-width indices to 
estimate stellar metallicities (Frogel et al. 200 \\, an d detailed 
abundance analysis of red giants (Orig lia et al.|2002"{|Origha~&| 
|Rich||2004] l, and supergiants ( |Davies et al. 2009a b ) at greater 
resolving powers in the //-band. 
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Fig. 5. Continuum signal-to-noise (S/N) results for simulated J- 
band spectroscopy (f exp = 10 hrs, R = 4000) from Table|4] 
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Fig. 6. Continuum signal-to-noise (S/N) results for simulated J- 
band spectroscopy (f exp = 10 hrs, R = 10000) from Table [4] 
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band spectroscopy (f exp = 10 hrs, R = 10000) from Table [5] 
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Fig. 10. Mean metallicities (and standard deviations) recovered from analysis of the simulated /-band spectra. Results are shown 
for the log[Z] =0.0 and -1.0 input templates, at R = 4000 and 10000, for both seeing values (0765 and 079) and configurations of 
natural guide stars (NGS; 'good' and 'poor'). 
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Fig. 9. As for Figure [s] but now for the template absorption-line 
spectrum, with Gaussian noised added (such that S/N = 50). 
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